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ABSTRACT 

To study the local Hubble flow, wc have run constrained dark matter (DM) simula- 
tions of the Local Group (LG) in the concordance ACDM and OCDM cosmologies, 
with identical cosmological parameters apart from the A term. The simulations were 
performed within a computational box of 64h~^Mpc centred on the LG. The initial 
conditions were constrained by the observed peculiar velocities of galaxies and posi- 
tions of X-ray nearby clusters of galaxies. The simulations faithfully reproduce the 
nearby large scale structure, and in particular the Local Supercluster and the Virgo 
cluster. LG-like objects have been selected from the DM halos so as to closely resemble 
the dynamical properties of the LG. Both the ACDM and OCDM simulations show 
very similar local Hubble flow around the LG-like objects. It follows that, contrary to 
recent statements, the dark energy (DE) does not manifest itself in the local dynamics. 
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.1^ ■ 1 INTRODUCTION 

5_] 1 It has been recently stated that the cosmological con- 
d ' stant (A) or its generalisation the dark energy (DE), 
manifests itself in the dynamics of the local universe 
(iBarvshev et al1l200l| . [Ch"ernin et al.ll2004 iTeerikorpi et al.] 



by means of constrained simulations (CSs, Kravtsov et al.l 
2C)0llKlvpin et all [200^ . [Martinez- Vagueroet al.ll2007l aiid 



Hoffman et al. 



2007l ~ 



20051 IChernin et all |2006| . IChernin et al.l l2007bl . aiid 



Chernin et al.ll2007d ). In these papers, the coldness of the 



local Hubble flow around the Local Group (LG) has been 
att ributed to th e exist ence DE. This has been supported 
by iMaccio et al.l ((20051) who analysed a set of N-body sim- 
ulations and concluded that indeed ...[their] results provide 
new, independent evidence for the presence of dark energy on 
scales of a few megaparsecs. These results, if correct, would 
have provided an independent corroboration to the DE com- 
ponent whose existence is otherwise inferred from observa- 
tions of distant objects and the early Universe. These au- 
thors used the term 'local' as describing the neighborhood 
of the LG out to a distance of a few Mpc. 

Much of the dynamical implications of the cosmolog- 
ical co nstant for the larg e scale structure were worked 
out by iLahav et al.l (|l99ll ). The LG constitutes a quasi- 
linear object and therefore its dynamics cannot be mod- 
eled by the linear theory or the spherical top-hat model. 
Consequently, we have recently studied the local universe 



The unique feature of the CSs is that 
their initial conditions are generat ed as constrained realiza - 
tions of Gaussian random fields (|Hoffman fc Ribaklll99ll ). 
The initial conditions are constrained by observational data 
and hence they are designed to reproduce the main gross 
features of the local large scale structure. As such they pro- 
vide the optimal tool for studying the dynamics of the LG, 
being a given individual but not an atypical object. In par- 
ticular the recent constrained flat A dominated ( ACDM) 
and open (OCDM) cold dark m atter N-body simulations of 
[Martinez- Vaguero et al.l ([2003) were designed to study the 
local dynamics in cold dark matter cosmologies with and 
without a DE component. These simulations are to be used 
here as a laboratory for thesting the hypothesis that the cold 
local Hubble flow is a signature of dark energy. We are less 
interested here in the actual coldness of the flow and more in 
the possibility that the DE affects the local flow. A thorough 
analysis of the issue of the coldness of the local flow is to be 
given elsewhere (Martinez- Vaquero et al, in preparation). In 
what follows 'local' is defined as the region contained in a 
sphere of radius R — 3 Mpc centred on the LG. 

The structure of the p aper is as follows. A very brie f 
review of the simulations of [Martinez- Vaouero et al.l (|2007l ) 
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is presented in Section [2] The selection criteria for LG can- 
didates are summarised in Section |3] The flow fields around 
the simulated LG candidates in the ACDM and OCDM sim- 
ulations are presented in Section |4l In Section [5]we compare 
the gravitational field around the LG candidates. A general 
discussion concludes the paper (Section [6]). 



CONSTRAINED SIMULATIONS OF THE 
LOCAL UNIVERSE 



Our C Ss have already been used in lMartinez-Vaauero et alJ 
(120071 ) and they are briefly summarised here. These are dark 
matter (DM) only simulations employing a periodic cubic 
computational box of 64h~^Mpc on a side using 256'^ par- 
ticles. Both models use the dimensionless Hubble constant 
of /i = 0.7 (where h = _ffo/100 km/s/Mpc), the power spec- 
trum normalisation erg = 0.9 and the cosmological matter 
density of fim = 0.3. The ACDM model corresponds to 
a flat universe with Qa = 1 — Q,m while for the OCDM 
model Q.A = 0. These cosmological parameters correspond 
to the so-called Concordance Model . We used the parallel 
TREEPM N-body code GADGET2 (|Springell , [2OO5I ) to run 
these simulations. For the PM part of the algorithm, we 
used a uniform grid of 512'' mesh points to estimate the 
long-range gravitational force by means of FFT techniques. 
The gravitational smoothing used to compute the short-scale 
gravitational forces correspond to an equivalent Plummer 
smoothing parameter of e = 15/i~^ kpc comoving. 

The number of particles used in these simulations (256^) 
provides a very mild mass resolution (1.3 x lO^h~^M0 per 
particle) which corresponds to a minimal mass of the DM 
halos of ~ 2.5 x 10^'^ h^^M©, for objects resolved with more 
than 20 dark matter particles. At such a resolution the inner 
structure of the main halos of the LG-like objects cannot be 
resolved, nor can the observed mass distribution of the LG 
nearby dwarfs be reconstructed. Yet, the dynamics on the 
scale of a very few Mpc is very well resolved. 

We set up initial conditions for these simulations in such 
a way that we can zoom in to any particular object with 
much more resolution. Thus, we generate the random re- 
alizations of the density fluctuation field for a much larger 
number of particles (up to 4096^). Then, we substitute the 
fourier modes corresponding to the small wavenumber by 
those coming from the constrained 256'' density field and 
make the displacement fields according to the Zeldovich 
approximation. Thus, we can now resimulate any particu- 
lar zone of the simulated volume with particles of variable 
masses, down to 4096 times smaller than the particle mass 
of the simulations used in this work. A comparison of the 
results of ACDM 256'' simulation with that from the LG- 
like systems resimulated at 4096^ resolution does not yield 
any significant differences in their Hubble diagrams (to be 
published). 

The algor ithm of constrained real izations of Gaussian 
random fields (|Hoffman fc Ribaklll99ll ) has been used to set 
up the initial conditions. The data used to constrain the ini- 
tial conditions of the simulations is made of two kinds. The 
first data set is mad e of radial velocitie s of ga laxies drawn 
from the MA RK HI (IWillick et aLlll997l ), SBF (|Tonrv et all 
I2OOII ) and the iKarachentsevI (2005|) catalogues. Peculiar ve- 
locities are less affected by non-linear effects and are used 



as constraints as if they were linear quantities (Zaroubi, 
Hoffman fc Dekel 199 9). T his follows the CSs p erformed by 
iKravtsov e'Ta^ (|2002l ) and iKlvpin et al.1 (120031 ). The other 
constraints are obtained from the catalog of nearby X-ra y 
selected clusters of galaxies (JReiprich fc Bohringeij I2OO2I ) . 
Given the virial parameters of a cluster and assuming the 
spherical top-hat model one can derive the linear overdensity 
of the cluster. The large scale structure, i. e. scales somewhat 
larger than 5h~^Mpc, of the resulting density and veloc- 
ity fields are strongly constrained by the imposed data. In 
particular all the resulting CSs are dominated by a Local 
Supercluster (LSC) - like object with a Virgo size DM halo 
at its center. The LG is not directly imposed on the initial 
conditions, but having reconstructed the actual large scale 
structure of the local universe a LG-like structure is very 
likely to emerge in the right place with dynamical proper- 
ties similar to the actual ones. The two simulations used 
here are based on the same random realization of the initial 
conditions. 



3 SELECTION OF LG-LIKE CANDIDATES 

The selection of LG candida tes is described in detail in 
[Martinez- Vaguero et al.l (120071). The selection of the objects 
is based on the iMaccio et al.l (|2005l ) criteria, which consist 
of: 



i. The group contains two MW and M31 like DM halos 
with maximum circular velocity in the range of 125 ^ Vc ^ 
270 km/s. 

ii. The two major DM halos are separated by no more 
than lh"^Mpc. 

iii. The relative radial velocity of the two main halos is 
negative. 

iv. There are no objects with maximum circular velocity 
higher than MW and M31 candidates within a distance of 
3h~^Mpc. 

V. The group resides within a distance of 5 to 12h^^Mpc 
from a Virgo like halo of 500 < 14 < 1500 km/s. 

DM halos a re found using both the Bound Density Max- 
ima alg orithm (iKlypin et"aU Il999l ) and the AMIGA Halo 
Finder (|GiU et al.l . l2004h . In the ACDM simulation 26 LG- 
like objects have been found and 43 in the OCDM one. Given 
the fact that both simulation are based on the same realiza- 
tion of the random Gaussian field we have identified 9 LG- 
like objects that appear in both simulations at about the 
same position and are very similar dynamically. We refer to 
these as the 'same' objects appearing in both simulations. 
These 'same' objects do not form any class by themselves 
and are statistically indistinguishable from the other LG- 
like objects. These objects are used here to exemplify the 
effect of the A term on the dynamics of the LG, as they are 
the same object evolving in two identical cosmologies and 
environments that differ only by their A term. 

The fact that there is no one-to-one coincidence of the 
LG-like objects of the two simulations should not be sur- 
prising. There are two reasons for that. First, the two cos- 
mologies are not identical and they differ in the linear grav- 
itational growth function. Second, the LG is a system in the 
quasi-linear regime and is far from being in dynamical equi- 
librium. Had we observed it at a slightly different time it 
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might not be qualified as a LG-like object according to the 
selection criteria assumed here. This is certainly the case for 
our simulated objects. Just a small miss match in the dy- 
namical phase of the objects between the two simulations 
can rule out an object in one or the other simulation from 
being a LG-like systems. 

In the present paper we are interested in comparing the 
local Hubble flow around LG-like objects. Providing that the 
selected systems fulfil all requirements their exact location 
is not important for the purpose of the analysis.. Therefore, 
we have used all the objects found within the computational 
box, regardless of their position with respect to the LSC. 
Some simulated LG-like objects reside close to the actual 
position of the LG but they seem to be dynamically indis- 
tinguishable from the others. 



4 THE LOCAL HUBBLE FLOW 

The local Hubble fiow around LG-like objects is probed by 
means of Hubble diagrams showing the radial velocities rel- 
ative to the objects center of mass within a distance of 3 
Mpc. In Figs. [T] and [2] we present the Hubble diagrams of 
4 randomly chosen candidates out of the 9 LG-like objects 
which appear in both simulations (Fig. [T] ACDM, Fig. [S] 
OCDM). The figures present all the DM halos around the 
chosen LG-like objects out to a distance of 3 Mpc. A careful 
comparison of the plots reveals that the Hubble diagrams of 
a given ACDM and OCDM simulated LG are very similar. 
In particular the r.m.s. value of the scatter around a pure 
Hubble fiow {an, assuming the true value of the Hubble con- 
stant of the simulation) does not vary statistically between 
the ACDM and OCDM cases. For the 4 objects shown in 
Figs 1 and 2, we find gh = 35, 38, 42 and 53 km/s for the 
ACDM objects and 41, 42, 55 and 59 km/s in the OCDM 
case. The values of an for the other 5 common candidates 
are 31, 51, 11, 41 and 79 km/s for the objects in the ACDM 
simulation and 61, 38, 18, 63 and 48 km/s in the OCDM 
one. 

Much of the theoretical expectations for the possible 
manifestation of the DE in the local fiow is based on the 
model proposed by IChernin et al.l l|2007d . and references 
therein). The model essentially assumes that the local grav- 
ity field around the LG can be decomposed into the con- 
tribution of the LG, modeled as a point particle, and the 
contribution of the DE: 



r (Mpc) Obs. ACDM OCDM 



, , GMlG , ^ rr2 

gppyr) = h ilf^Hor 



(1) 



The zero gravity surface is defined by gpp{Rv) = 0. The 
radius of the zero gravity surface, Rv, plays a critical role 
in that simple model. A central prediction of the model is 
that the local Hubble fiow should not contain galaxies with 
radial velocities smaller than the escape velocity (see the 
Appendix), calculated under the assumption that the grav- 
itational field is given by the point particle approximation 
(jChernin et al.ll2007a ). This prediction excludes galaxies re- 
siding within the LG itself, namely within 0.7 Mpc. To test 
the prediction the radial escape velocity profiles (Eq.[6} have 
been plotted in both Figs. [T] and [2] as solid (ACDM model) 
and dashed (OCDM model) lines. From here on the term 
'escape velocity' refers to the one calculated under the as- 
sumption of the point particle approximation. 
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Table 1. The value of an (in units of km/s) of the LG, compiled 
from the Karachentsev data, and of the ACDM and OCDM LG- 
hke objects combined together, in the manner of Fig. |3] Two 
distance cuts are used for calculating an- 



Inspection of the ACDM Hubble diagram (Fig.[Tl ) shows 
that indeed the prediction of lChernin et al.l (|2007bl ) is con- 
firmed: only two LG-like groups have, within the range 
(0.7 — 3) Mpc, a very few halos each with a peculiar velocity 
smaller than the ACDM escape velocity. However, this be- 
haviour is reproduced by the OCDM LG-like objects equally 
weU (Fig.IH). 

To increase the statistical significance of the Hubble di- 
agram analysis we have considered all the LG-like objects 
in the ACDM and OCDM simulations. This is performed 
by plotting the radial velocities of all halos near the LG-like 
groups against their distance r. Again, the Hubble diagram 
of all LG-like objects in both cosmologies looks very simi- 
lar. In fact, the fraction of halos below the escape velocity is 
somewhat smaller in the OCDM objects than in the ACDM 
ones. We conclude that the ACDM escape velocity predic- 
tion is reproduced by the OCDM simulation. 

The paper focuses mainly on the possible role of the 
DE in the dynamics of the LG. A thorough analysis of the 
coldness of the local flow will be given elsewhere (Martinez- 
Vaquero et al, in preparation). Here a very brief summary 
of the subject is gi ven. The very local Hubble flow has been 
recently studied bv lKarachentsev et al.l l)2007h and their cur- 
rently updated catalog of local peculiar velocities has been 
analyzed here (I. Karachentsev, private communication). 
Table I presents the value of gh taken over all the DM ha- 
los (simulations) or galaxies (data) in the range of [0.75 — 2] 
Mpc and [0.75 — 3] Mpc of the Karachentsev's data and of 
the ACDM and the OCDM LG-like objects. The cumulative 
distribution of an (calculated over the range [0.75 — 3] Mpc) 
is presented in Fig. |4] The plot shows that more than half 
the LG-like objects in both models have a an ^ 60 km/s. 
So, many objects have a fiow as cold, or colder, a s the ac- 
tual LG. Yet, as was pointed bv lMaccio et al.l ((20051) the real 
problem of the coldness lies with the relation between an 
and the mean density around the objects. 

It follows that the local Hubble fiow around ACDM 
and OCDM LG-like objects is essentially indistinguish- 
able. Th is stands in c l ear c o ntradiction wi t h pre vious 
claims oflBaryshev et"aLl (|200lf ). [Chernin et al.l (|2007cl ') and 
IChernin et al.l l|2007ah . Also, both the ACDM and OCDM 
LG-like groups obey equally well the escape velocity predic- 
tion of the fiat -A cosmology, as if they are not affected by 
the A term. 



5 THE LOCAL GRAVITATIONAL FIELD 

To understand the possible reason for the discrepancies be- 
tween the present results and the and the model predictions 
we have studied the nature of the local gravitational field. 
The prime motivation here is to check the validity of the 
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IChernin et alj (|2007d ) model of the gravitational field (Eq. 
[T]), which corresponds to the full gravitational field expressed 
in physical, and not co-moving, coordinates. The relation 
between the peculiar gravit y (output of GADGET) a nd th e 
physical one was derived bv lMartinez-Vaguero et alJ ((20071), 
but it is repeated here for the sake of completeness. 

The physical r and comoving x coordinates are related 
by: r = ax. The gravitational field equals the physical ac- 
celeration of an object r = g. The GADGET code provides 
an acceleration-like term defined as: 



Id 
a at 



Vp + -H" ■ Vp 



(2) 



where Vp = r — Hr is the peculiar velocity, H is Bubble's 
constant and a is the expansion scale factor. It follows that 



i^ = fp -I- r- = fp -f 



]^^m+^a]H'"v. 



(3) 



Namely, the linear term corresponds to the unperturbed 
Friedman solution and fp to the fluctuating componenllj. 

The gravitational field is taken with respect to the LG 
reference frame. So that, one finally obtain: 



5 = (fp-fp^'')i + 



--^A/ + f^A 



H' 



(4) 



where r is the distance from the center of mass of the LG. 
This is the field acting on each dark matter particle. The 
total acceleration of halos was computed by averaging this 
quantity over all particles belonging to each halo. The accel- 
eration is scaled by Hq x 1 Mpc. In such scaling the unper- 
turbed gravitational acceleration of a shell of radius 1 Mpc 
equals to —go, where go is the cosmological deceleration pa- 
rameter. 

In Figs. [5] and |6] we compare the radial component of 
the exact (i.e . in the sense of the si mulations) gravitational 
field with the lChernin et al.l (|2007d ) model (gpp), as traced 
by the DM halos around the LG-like objects. Only the fiuc- 
tuating component of the gravitational field is shown in the 
figures, namely for the numerically exact field it is the radial 
component of fp and for the point particle approximation it 
is fp^pp = -GMLc/r^ + ^MH^r/2. We show aU DM ha- 
los in the distance range of 0.7 ^ r ^ 3.0 Mpc from the 
same LG-like objects that were presented in Figs. [T] and [21 
Since DM halos closer that r = 0.7 Mpc are affected by the 
two-body dynamics of the LG, they are excluded here. Fig. 
[7] shows the scatter plot of the gravitation field of all the 
LG-like objects of the ACDM and OCDM simulations. The 
plots show that the numerically exact value of the radial 
component of fp and the point particle prediction {fp,pp) 
are very poorly correlated. A linear regression analysis finds 
a correlation coefficient of 0.40 (0.23) and a slope of 0.35 
(0.16) for the ACDM (OCDM) model. It is clear that the 
point particle model fails to reproduce the actual gravity 
field, and therefore it cannot account for the dynamics of 
the LG. 



1 Note the typo in Eq. A3 of [Martinez- Vaguero et al.l bOO?! ') 
where Ha was omitted. However, this did not affect the result, as 
only the fluctuating term of g was considered there. 



6 DISCUSSION 

The most striking result of this paper is that the local Hub- 
ble fiow around LG-like objects in the OCDM model is al- 
most indistinguishable from the ACDM fiow. To the extent 
that the models do differ it is the OCDM model that has 
somewhat colder Hubble flow than the ACDM one. It fol- 
lows that the local flow is not affected by the DE and does 
not manifest the present epoch dominance of the DE. 

One should not be surprised by the departure of the sim- 
ulated local Hubbl e fiow from the predic tion of the simple 
model proposed bv lChernin et al.l (|2007cr ). First, the actual 
gravitational field deviates considerably from the predicted 
one. Second, the gravitational dynamics is not local and the 
tidal field plays an import ant ro l e in the quasi-linear regime 
(Hoffman 1986. Hoffmanl | l989| . IZaroubi fc HoffmanI Il993l . 
Ivan de Weveaert fc Babullll994l . [Del Popolo et al.ll200ll )~It 
follows that the dynamics depends not only on the local field, 
but is also affected by the shear, namely the tidal field. The 
shear breaks the simple linear relation of the density and 
velocity fields of the linear regime and therefore the local 
density field cannot account for the local Hubble fiow. 

The comparison of the ACDM and OCDM simulations 
shows that they yield very similar LG-like objects with virtu- 
ally identical local Hubble fiows. It follows that the dynam- 
ical properties of LG-like objects and their environments, in 
the linear and quasi-linear regime, depend mostly on the cold 
matter content of the universe, namely ilm, and only weakly 
on the DE. This is another manifestation of the fact that 
the properties of the cosmic web, expressed in co-moving 
coordinates, depend mostly on fJm and hardly on the DE 
([Hoffman et al.|]2007l ). 
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APPENDIX 

Assuming the local gravitational field around the LG is in- 
deed given by the point particle approximation (Eq. [TJ , the 
effective Newtonian potential is given by 



GMlg ftAHy 



(5) 
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The effective Newtonian energy (per unit mass) is simply 
given as e = tj^/2 + 0(r). In the presence of the A term the 
potential reaches a maximum at Rv at which the potential 
peaks at ev = 4>[Rv)- It follows that a particle is unbound 
to the LG if its energy is larger than ev- The escape velocity 
is therefore given by 



In the case of a vanishing cosmological constant the classi- 
cal expression of Eq. [5]and|6]is recovered upon substituting 
Rv = oo and SIa = 0. It should be reemphasised here that 
the present derivation is done under the assumption of the 
point particle approximation. The analysis of the simula- 
tion proves that the assumption is inapplicable to the LG 
systems. 
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Figure 1. The Hubble diagrams around four LG-like objects in the ACDM simulation. The scatter plots represent the radial peculiar 
velocity of the DM halos vs. the distance from the MW and M31-like DM halos. The solid curve corresponds to the escape velocity 
profile of the ACDM model, calculated under the assumption of the point particle approximation. For reference the escape velocity of 
the corresponding OCDM model is given as well (dashed line), namely it is calculated as if the A term is missing. The value of au, taken 
over the range 0.75 ^ r ^ 3 Mpc and Ry of each object is given. 
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Figure 2. The Hubble diagrams around four LG-liltc objects in tlie OCDM simulation. The four LG-like objects shown here are the 
OCDM counterparts of the ACDM ones shown in Fig. [T] The structure and notations of the plots are the same as in Fig. [T] 
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Figure 3. Combined Hubble diagram of 26 LG candidates in the ACDM (left panel) and 43 candidates in the OCDM (right panel) 
models. The radial distance r is scaled by the value of Ry of each object. The values of au within 2 and 3 Mpc are given in Table 1. 
The escape velocity curves are plotted in the same manner as in Figs. [T] and [2] < Ry > is the mean Ry of all the LG-like objects for 
each simulation. 
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Figure 4. The fractional cumulative uh function of LG-like objects, namely the fraction of objects with uji lower than a certain value. 
Full line corresponds to the 26 ACDM objects and the dashed one to the 43 OCDM objects. The dispersion around a pure Hubble flow, 
aii^ is calculated over the range 0.75 SC r sC 3 Mpc. 
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Figure 5. The ACDM gravitational field: A scatter plot of the exact gravitational field g (Eq. [4]| vs. the approximated one gpp (Eq. 
[TJ experienced by DM halos around LG-like objects. The four objects and panels correspond to the ones in Fig.[T] Only DM halos in the 
range of 0.75 ^ r ^ 3.0 Mpc are plotted here. The gravitational field is scaled by Hq x 1 Mpc. 
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Figure 6. Same as Fig. [5]but for the OCDM simulation. The LG-like objects and panels correspond to the ones in Fig. [2] 
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Figure 7. Scatter plot of the fluctuating component of the gravitational field predicted by the point particle model {fp^pp) against the 
exact value calculated by the simulations. The left panel shows the distribution of the DM halos in the vicinity (0.75 ^ r ^ 3.0 Mpc) of 
the 26 LG-like ACDM objects and the right panel exhibits the 43 OCDM objects. The gravitational field is scaled by Hq x 1 Mpc. 



